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Simulation of high-intense beam transport in electrostatic accelerating column* 
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An electrostatic accelerating column was designed and fabricated by Lanzhou University for an intense D- 
T/D-D neutron generator. In order to achieve a neutron yield of 5.0 x 10'* n/s, a deuteron beam of 30 mA, 
accelerated to 400 kV, and transported in the electrostatic accelerating column smoothly are required. One 
particle-in-cell code BEAMPATH was used to simulate the beam transport, and the IONB1.0 code was used 
to simulate the intense beam envelopes. Emittance growths due to space charge effect and spherical aberration 
were analyzed. The simulation results show that the accelerating column can transport deuteron beam of 30 mA 
smoothly and the requirement for the neutron generator is satisfied. 
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I. INTRODUCTION 


In recent years, development of high-current ion accelera- 
tor have evoked considerable interest for various applications 
like accelerator-driven systems (ADS) [1], accelerator-based 
neutron sources [2], next-generation radioactive ion beam fa- 
cilities [3], ion cancer therapy facilities [4, 5], etc. A D- 
T/D-D neutron generator which is a monoenergetic neutron 
source driven by an electrostatic accelerator is being devel- 
oped by Lanzhou University. An electrostatic accelerating 
column for accelerating and transporting high-intense beam 
was designed and installed on the neutron generator. In or- 
der to achieve a neutron yield of 5.0 x 10 n/s, a deuteron 
beam of 30 mA, accelerated to 400 kV and transported in the 
electrostatic accelerating column smoothly, is required [6, 7]. 


This paper aims to analyze and simulate the transport of 
high-intense deuteron beam in the electrostatic accelerating 
column for the D-T/D-D neutron generator. BEAMPATH [8], 
a well proven code using the PIC method, is used to carry out 
the beam transport simulations and investigate the emittance 
growth due to space charge effect [9] and spherical aberra- 
tions [10] in the electrostatic accelerating column, and a com- 
puter code named IONB1.0 for simulating intense beam en- 
velopes which based on K-V equation has been developed 
and used. 
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I. ANALYSIS OF BEAM EMITTANCE GROWTH 
A. The accelerating column 


The accelerating column is 1370 mm in length (Fig. 1). It 
consists of the accelerating zones of Al and A2, both being 
the size of P90 mm x 110mm. They can be seen as two elec- 
trostatic lenses. Two electrodes biased negatively are installed 
on the accelerating column exit, for repelling secondary elec- 
trons back into the accelerating column and keeping them free 
from the ion source. Between the two electrodes is a per- 
manent magnet lens at ground potential. The magnet lens 
and electrodes form a space charge lens [11], so as to pre- 
vent beam diverging excessively, and adjust the size of beam 
envelope at the accelerating column exit. 
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Fig. 1. (Color online) Schematic diagram of the accelerating colum- 
n. 


B. Effects of spherical aberrations 


In general, transport of space charge dominated beam is 
analyzed in linear approximation method. In the presence of 
linear focusing field, the beam emittance remains constant. 
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Realistic focusing elements possess strong aberrations, result- 
ing in distortion of phase space area occupied by the beam. 
Among others, the spherical aberration cannot be eliminated, 
hence the most significant effect on particle dynamics. The 
spherical aberration coefficient C; is defined as [10, 12] 


C, = (f/r)? Ar, (1) 


where, f is the local length, r is the radial distance of an 
incident ray parallel to the axis, and Ar is the radial distance 
at the paraxial-image plane. 

The change of slope of particle trajectory can be written 
via focal length f and aberration coefficient C, [10, 12]. 


A(dx/dy) = —(x/f) [1+ Cs(r/f)?/f] - (2) 


Assume the initial beam distribution is described by the 
ellipse 
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where, R is the envelope size of beam, and € is the beam 
emittance. 

The transformation from initial particle variables before 
lens (xo, xo) to that after the lens (x, x’) is given by 


£ = £0, (4) 
Q 
z = xh (1+ rs) - (5) 
Change variables (x, x’) to action-angle variables (J,W) 
Ave = V2J cos, (6) 


R 
(« F z) * -= aT sinw. (7) 
fJ Ve 
Then beam ellipse distortion: 
T + T?2w sin Y cos? Y + T?v? cos? Y = 1, (8) 


where, T = 2J/e and v = C,R*/(ef*). Without nonlinear 
perturbation, v = 0, the shape of beam emittance is round. 
While phase space areas occupied by beam before and after 
lens are the same, the effective area, occupied by the beam, 
is increased. The beam emittance can be estimated as a total 
area of elements dada’ occupied by the beam. Actual areas 
of the beam in both cases are the same, but the number of 
elements covered by the beam at phase pane is different. The 
increase of effective beam emittance as a square of product of 
minimum and maximum values of T: 


Eeff/E = (Taastan) (9) 


where Tmax and Tmin are determined numerically. Depen- 
dence of emittance growth versus parameter v can be approx- 
imated by 


ear/e = (1+ Ko), (10) 
where K ~ 0.4. Finally, effective beam emittance growth 
due to spherical aberrations can be approximated by 


Eeff/E = E +K (aRt ae (11) 
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C. Aberrations of space charge effect 


Space charge density and space charge field of the beam of 
Gaussian distribution in free space are given by [13, 14] 


21 r2 
= g0 12 
p(ro) Eie exp ( ia). (12) 


hel 910 (13) 
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where, €o is permittivity of vacuum; J is the beam current; 
and 3 = u/c, vis the particle velocity and c is the light speed 
in vacuum. 

Nonlinear function in space charge field can be expanded 
as 


re re rå 
f(ro) =1 exp ( 2 )- 2— +. (14) 
Ro Ro 
At the initial stage of beam emittance growth we can as- 
sume, that particle radius does not change, while the slope of 
the trajectory changes. Then, we have the following nonlinear 
transformation: 


r = To, (15) 


Fai 22P? E 2P? a 


RB “Re TO (16) 
where, P? = 21/(I.837?) is the generalized perveance, ~y is 
relativistic factor, and I, = 47tegmc? /q is the characteristic 
beam current. 

Use the same way as above, effect of spherical aberration 
on beam emittance can be determined as 


4 Iz 
v= aa 
Bee Iae 


(17) 


Therefore, space charge induced beam emittance growth in 
free space is 


cent /E = [1 +s [I/ (LBN (z/e)"| k , (18) 


where, z is the drift length and x = 0.6 is beam uniformity 
for Gaussian beam. 


I. ENVELOP SIMULATION BASED ON K-V EQUATION 
A. High-current ion beam envelope equation 


For designing the intense neutron generator, a comput- 
er code of IONB1.0 was developed. Assuming a continu- 
ous, elliptically symmetric particle beam, its dynamics can be 
modeled by the K — V (Kapchinsky—Vladimirsky) coupled- 
envelope equations [13, 15, 16]. Let the z coordinate rep- 
resents the position along the design trajectory, and the x-y 
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plane is the transverse plane for the particle beam, we have 
the high-current beam envelope equation [17]: 


/ n oT c2 Yo 
x'a E a XIEX zor 
T 2p 4y = X+Y Xp 
1 n" DANI e2 Yo A 
Y” +2 yY' 4+ Y +K,Y w’ o 
= 2p T 4y Tan X+Y Y’ọ 


(19) 
where, X (z) and Y (z) are dimension of the semi-axes of 
beam envelope in the x and y planes, respectively. The prime 
indicates differentiation with respect to z. y, y’ and wy” are 
the generalized potential, the electric field strength and the 
derivative of electric field strength at z axis (Beam axis), re- 
spectively. Ką, and K, refer to the focusing functions of 
transport elements in the x and y directions respectively. yo 
is the normalized potential, €z, and €,, are the normalized 
emittances of beam on the (x, x’) and (y, y’) phase planes re- 
spectively. 2II is the generalized perveance, it reflects space 
charge effect of the high-intense beam, and can be calculated 
by following formulas 


=" (2n)- 2973, (20) 


where, J is the beam current, F is the factor of space charge 
compensation, 77 is the charge-to-mass ratio of charged parti- 
cle, and £ọ is the permittivity of vacuum. 


B. Numerical simulation method 


For calculating X, the K — V equation can be expressed as 


2JI €2 E0 
+ 
X+Y X’ọ 


/ {i 
X" =- a 


. (21 
2 Ee (21) 


As in Fig. 2, the z axis is divided into a series of micro- 
units, Az. When Az is tiny enough, we may approximately 
hold that X” is a constant A in each micro-unit. 


Fig. 2. Divided diagram for the beam axis. 


After the beam passes through a Az, its envelope and di- 
vergence may be approximately calculated by Eq. (22) 


X'=AAZ+B 


Donet i (22) 
X= z447 + BAZ +C 


where, A, B and C can be fixed by the initial conditions of 
beam envelope. As in Fig. 2, if the initial conditions at point 
P are (Xp, Xb), in that case, B = Xb and C = Xp. A 
can be calculated by Eq. (21). The beam parameters at next 
interval point Q(XQ, Xo) are determined by Eq. (22). Then, 
Q is regarded as the start point, we can calculate the beam 
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parameters at next point R( Xp, Xp). So long as Az is tiny 
enough, the better consults of beam parameters can be got 
at each interval point. The same method is used to calculate 
the Y. 


IV. BEAM SIMULATION 
A. Optimization accelerating column fields 


The distribution of electric fields is not only to meet the re- 
quirements of increasing beam energy and inhibiting beam di- 
vergence excessively, but also to make 400 kV voltage have a 
reasonable allocation to avoid local field strength be too large 
to produce discharge. Using the POISSION/SUPERFISH 
code [18, 19], the distribution of equipotential surface was 
simulated as shown in Fig. 3, the electric field strengths in all 
regions are less than half of its maximum value of the break- 
down electric field strengths, the electric field directions are 
marked by the arrows. Figure 4 shows the map of radial fields 
E,(z,r) and E,(z,r) at 0—4 cm radii. 


Fig. 3. (Color online) Equipotential surface distribution of acceler- 
ating column. 


B. Initial parameters 


There are a number of PIC codes developed for simulat- 
ing space charge-dominated beam dynamics [20, 21]. Among 
them BEAMPATH can simulate unbunched continuous beam 
dynamics using the PIC method. Our simulation began with 
the point S (Fig. 1). A Gaussian distribution was assumed 
for a 20keV axial symmetric deuteron beam. According to 
experimental data of the duoplasmatron ion source [22], ini- 
tial parameters at S of deuteron beam were chosen as follows: 
beam envelope radius, X = Y = 10mm; divergence angle, 
X’ = Y’ = 50 mrad; natural emittance, € = 150mm mrad; 
and normalized beam emittance, 0.7 mm mrad. The E(z,r) 
field map of the accelerating column was calculated using 
the POISSON/SUPERFISH, and inputting into the BEAM- 
PATH and IONB1.0, respectively. Beam intensity was set 
to 30mA and 0mA for considering and ignoring the space 
charge effect, respectively, when BEAMPATH was used. For 
IONB1.0, the space charge compensation had a significant 
impact on simulation results. In our previous study, fully 
unneutralized beam was assumed in the accelerating zones, 
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Fig. 4. Radial field maps of (a) Æ, (z, r) and (b) Ez (z, r). 


with the space charge compensation factor of 0.8 in the non- 
accelerating zones [17], but it should be 0.65 in the non- 
accelerating zones and 0 in the accelerating zones, compar- 
ing the simulation results of IONB1.0 with the BEAMPATH 
results. 


C. Simulation results 


The charged particle motion in electrostatic accelerating 
column depends on the external electrostatic field and the 
space charge self-field, and the space charge force causes 
beam divergence and emittance growth. The simulation re- 
sults indicate that the deuteron beam of 30 mA can go through 
the accelerating column successfully. 

The simulation results with space charge effect are shown 
in Fig. 5(a). The maximum of beam envelope radius at accel- 
erating column exit is about 2.5cm using IONB1.0, and the 
beam envelop radii of the vast majority of deuterons at ac- 
celerating column exit are 2-3 cm using BEAMPATH; while 
without space charge effect (Fig. 5(b)), the maximum of beam 
envelope radius at accelerating column exit is about 1.2cm, 
and the beam envelop radii of the vast majority of deuterons at 
accelerating column exit are —1.2cm. The beam divergence 
and the maximum of beam envelope radius are significantly 
less than those in Fig. 5(a). 

Figure 6(a) illustrates the oscillation of the emittance along 
the accelerating column of both deuteron beam of 0 mA and 
30mA. Significant beam emittance growth can be seen in 
the BEAMPATH simulation. The beam emittance of 0 mA 
differs from that of 30 mA in the accelerating zones, and they 
vary differently in different accelerating zones. The beam 
emittance at the accelerating column exit is 0.3cm mrad, be- 
ing 4.3 times greater than the original emittances. As shown 
in Fig. 6(b), the space charge effect plays a significant role 
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Fig. 5. (Color online) Envelopes along the accelerating column (a) 
with (30 mA) and (b) without (0 mA) space charge effect; the lines 
with red circle represent the simulation envelopes used IONB1.0, 
and the solid black lines represent the simulation envelopes used 
BEAMPATH. 


in the beam emittance growth in the accelerating zones, the 
maximum of €39/€o, which is the beam emittance ratio of 
30mA to OmA, is 1.24 in the accelerating zones, but it is 
nearly 1 in the non-accelerating zones. This shows that the 
beam emittance increased rapidly before the beam incident to 
A1, inhibited in A1, and decreased in A2. 
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Fig. 6. Variation of beam emittance along accelerating column (a) 
and beam emittance ratio of 30 mA to 0 mA (b). 
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(b) Phase space, 30 mA 
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Fig. 7. Deuteron distributions in phase space at the accelerating column exit, (a) without (0 mA) and (b) with (80 mA) space charge effect; 
and in real space at the accelerating column exit, (c) without (0 mA) and (b) with (30 mA) space charge effect. 


Figure 7 shows the final deuteron distributions in phase s- 
pace or real space, without or with the space charge effect. 
With the space charge effect, the beam divergence increas- 
es (Figs. 7(a) and 7(b)) and the hollow beam is formed 
(Figs. 7(c) and 7(d)). An annular distribution is formed by a 
large number of deuterons, with the inner radii of about 2 cm 
and the outer radii of about 3 cm (Fig. 7(d)) at the accelerating 
column exit, indicating that the electrostatic accelerating col- 
umn can transport space charge-dominated beam of 30 mA, 
as required by the design. 
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